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Abstract: We present the generating function approach to the perturbative exponentia¬ 
tion of correlators of a product of Wilson lines and loops. The exponentiated expression is 
presented in closed form as an algebraic function of correlators of known operators, which 
can be seen as a generating function for web diagrams. The expression is naturally split 
onto two parts: the exponentiation kernel, which accumulates all non-trivial information 
about web diagrams, and the defect of exponentiation, which reconstructs the matrix ex¬ 
ponent and is a function of the exponentiation kernel. The detailed comparison of the 
presented approach with existing approaches to exponentiation is presented as well. We 
also give examples of calculations within the generating function exponentiation, namely, 
we consider different configurations of light-like Wilson lines in the multi-gluon-exchange- 
webs (MGEW) approximation. Within this approximation the corresponding correlators 
can be calculated exactly at any order of perturbative expansion by only algebraic manip¬ 
ulations. The MGEW approximation shows violation of the dipole formula for infrared 
singularities at three-loop order. 
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1 Introduction 

The dynamics of Wilson lines and loops is responsible for nearly every aspect of gauge 
theories. Naturally, the correlators of Wilson lines and loops are one of the most attractive 
objects of theoretical investigations. In this article we present the generating functional 
approach to the exponentiation of a general correlator of Wilson lines. 

The exponentiation of a correlator of Wilson lines is the fundamental property of the 
perturbative approach. As the name suggests, the exponentiation property allows one 
to present the perturbative series for the correlator as an exponent. The non-triviality 
of this statement consists in the special diagrammatic of the exponentiated expression. 
Namely, only connected diagrams enter the argument of the exponent. Moreover, often 
the exponentiated diagrams possess some additional non-trivial properties, being a subset 
of the initial (not exponentiated) series of diagrams. Such an approach can be considered 
as an exact resummation of the initial diagram series. In this way, the knowledge of the 
rules for exponentiated diagrammatic significantly simplifies analysis of the Wilson lines 
correlators within the perturbation theory. 
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The exponentiation property for Abelian gauge theories has been understood a long 
time ago [1]. The exponentiation for a non-Abelian gauge theory is significantly more 
complicated, and until the recent time, only the exponentiation property for the cusp 
configuration of Wilson lines and for a Wilson loop was known [2-4] . The exponentiated 
diagrams for the cusp were called “webs” [2], and were investigated in great details. They 
play an important role in the description of the strong interaction dynamics. However, 
nowadays more complex configurations of Wilson lines are interesting. During the last 
decade several approaches to the exponentiation of a general correlator of Wilson lines 
were presented [5-8] (the brief review of these approaches is given in sec. 2). 

In this article we present and elaborate details of the exponentiation approach via the 
generating function, which has been partially presented in [9]. The main object of our 
interest is the correlator of any number of Wilson lines of any paths and of any group 
representations 

5iv({7i,/i},-,{7iv,/w}) = X ... X (1.1) 

where $ is a Wilson line of representation / along the path 7 : 

= Pexp dr 7 ''('r) 2 l“ ( 7 (T))tL^'^ , ( 1 . 2 ) 

where 7 (r) parameterizes the path, and 7 is a tangent to the path at point r, are 
the generators of the gauge group in the representation /. We stress that the Wilson line 
(1.2) is a matrix in the color space, and hence, the correlator (1.1) is a multi-matrix. The 
disentangling of the matrix structure is the main difficulty of the exponentiation in a non- 
Abelian gauge theory. In particular, inaccurate treatment of the matrix structure gives rise 
to the incomplete exponentiated expression in ref. [9] . This issue is corrected in the present 
article. 

The method of exponentiation presented in the article, which we call the generating 
function exponentiation, is novel and qualitatively differs from the methods presented in 
the literature. The main feature of the generating function exponentiation is that the 
argument of exponent is presented in closed form, namely, as a function of correlators 
of certain operators. This function has the meaning of the generating function for web 
diagrams. 

Technically, the main difference between the generating function exponentiation and 
other methods of exponentiation consists in the splitting of the problem onto two principally 
different mathematical tasks: the exponentiation of a scalar operator in a non-Abelian field 
theory, and the exponentiation of a matrix object. Separately these tasks are effortless. 
However, together they produce many various diagrams and factors, the straightforward 
exponentiation of which is a cumbersome work. Therefore, despite many parallels with 
exponentiation methods presented in literature, we found the final result unique. The 
detailed explicit comparison of the exponentiation approaches is given in sec. 5. 

The article is composed as follows. In sec. 2 we make a brief review of exponentiation 
methods. The main aim of this section is to introduce a minimal terminology. In sec. 3 we 
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present elementary introduction to the generating function approach to the exponentiation, 
in particular, we demonstrate the exponentiation in the Abelian gauge theory and the 
exponentiation of real exchanges. Section 4 is the main section of the article. In this 
section we present the generating function approach to the non-Abelian exponentiation 
and derive all the main formulae. Finally, sec. 5 presents a set of explicit applications of the 
generating function exponentiation. In particular, in sec. 5.1 on the example of the two- 
loop cusp, we perform diagram-by-diagram comparison of the presented approach with 
existing approaches, while in sec. 5.2 we demonstrate the effectiveness of the method by 
evaluating some configurations of light-like Wilson lines in the multi-gluon-exchange-webs 
approximation. 

2 Brief review of exponentiation methods 

In the paper we present an exponentiation method for a product of Wilson lines and 
compare it with other exponentiation approaches presented in literature. The methods 
(the one presented here and those taken from literature) have similar and distinct points. 
In order to clarify the presentation let us shortly describe the exponentiation methods and 
introduce necessary terminology. 

The most straightforward approach to the exponentiation is to consider the logarithm 
of the perturbative series. Performing the perturbative expansion of the logarithm and 
combining diagrams together one obtains the exponentiated perturbative series. Within 
the composition of diagrams many terms cancel and the remainder is called web diagrams, 
or webs [2]. We name such an approach to the exponentiation as the diagrammatic expo¬ 
nentiation. 

For the cusp and for a Wilson loop the diagrammatic approach gives the following 
result for the exponentiated series of diagrams. The color factor of the diagrams (of the 
non-exponentiated series) should be replaced by the modified color factors [3, 4], 

tr($) = ^ C{d)F{d) = exp ( ^ (7((i).F(d) J , (2.1) 

d \ d ) 

where C{d) and F{d) are the color factor and the kinematical part of a diagram d, respec¬ 
tively. The modified color factor C is obtained by the recursive procedure 

C{d) = C{d) - n (2.2) 

d' 

where the sum runs over all decompositions of d into two-Wilson-lines-irreducible subgraphs 
w (webs). In the consequence of the expression (2.2) only the two-Wilson-lines-irreducible 
diagrams have non-zero modified color factors. These diagrams form the set of webs for 
the cusp, or for the Wilson loop. 

For a general configuration of Wilson lines the diagrammatic exponentiation has been 
considered in ref. [7]. It has been shown that for the general case no simple selection 
criterium exists, and all diagrammatic topologies enter the exponentiated series. 
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Another method of exponentiation was suggested in ref. [5] and elaborated in [6, 8, 10]. 
The method is based on the replica trick and hence, we name it as the replica exponentia¬ 
tion. 

Briefly, the course of the replica exponentiation is the following. In the first step, 
the correlator of (replica-ordered) Wilson lines is considered in the replicated theory, i.e. 
the theory that consists of Wep copies of the original theory (replicas). The obtained 
diagrammatic expansion depends on Wep- In the second step, the terms proportional to 
Wep are taken, and the rest terms are discarded. These terms form the argument of the 
exponent in the original theory. The resulting expression can be presented in the following 
form 


(<h X <h... X <h)/j = exp 


Y,Hd)Rd,d'C{d!) 


d.d' 


IJ 


(2.3) 


where R is the web-mixing matrix [6, 10], and I and J are the color multi-indices. The 
product RC generalizes the concept of the modified color factor of the color-singlet case 
( 2 . 2 ). 

The replica exponentiation gives deeper theoretical understanding of the exponenti¬ 
ation procedure, and allows one to make certain general conclusions about webs. For 
example, using the replica exponentiation it has been shown that color factors of the 
exponentiated expression contain only the color connected combinations of diagrams [8]. 
Recently, the method has been improved by the effective vertex formalism [8], and in this 
form it was used for analysis of particular multi-loop webs [11]. 

The disadvantage of the replica method is that it does not present any universal ex¬ 
pression for the series of web diagrams. The web-mixing matrix should be calculated at 
every order of the perturbative expression independently, what makes any general con¬ 
sideration difficult. In this aspect, the replica exponentiation resembles the diagrammatic 
exponentiation for the product of Wilson lines [7] (for explicit comparison see the appendix 
of ref. [10]). 

In this article we present a novel exponentiation approach. The main distinctive feature 
of the approach is the presentation of the exponentiated series in closed form, namely, as a 
function of correlators of known operators. This function can be interpreted as a generating 
function for web diagrams [9]. Thus, we name the presented method as the exponentiation 
via generating function or, for shortness, GF exponentiation. 

Although GF exponentiation is obtained independently from the diagrammatic expo¬ 
nentiation and the replica exponentiation, these approaches have many common points. 
GF exponentiation can be viewed as a union of both approaches. As we show later, the 
generating function consists of two parts. One of those can be easily traced in the replica 
method (in its formulation via effective operators [8]), while another is closely related to 
the diagrammatic approach. However, in contrast to the diagrammatic and the replica 
exponentiations, the GF exponentiation has a visual and simple expression. 
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3 Generating function approach to exponentiation 

In this section we discuss the basics of GF exponentiation. The material of the section is 
somewhat trivial and generally can be collected from textbooks on quantum field theory. 
Therefore, the main aim of this section is not to present a new material, but to introduce 
notations, and the general way of approach. 

We also demonstrate GF exponentiation in Abelian gauge theory. This presentation 
serves as a pedagogical example, which is later generalized on the non-Abelian case. The 
presentation of this part closely follows the presentation in ref. [9]. 

In the last subsection, we briefly discuss the generalization of GF exponentiation on 
the correlators with real exchanges, which are very important in practice. 


3.1 The foundation of exponentiation 

In ref. [9] it was shown that the quantum average of any operator, which has the form 
of an exponent of another operator, can be presented as an exponent of the series of 
connected diagrams. In fact, this statement is a reformulation of the famous relation 
between the partition function and the series of connected diagrams (3.2). The relation 
holds in any quantum field theory, therefore, in this section we use maximally abstract 
notation, denoting all fields of a theory as A, and do not specify their quantum numbers. 
The partition function of a quantum theory reads 

Z[J] = j da Jix)o{x)^ (3 

where 0{x) is a composite operator, J{x) is the source, and ^[A] is the action of theory. 
The diagrammatic expansion of Z[J] consists of all possible Feynman diagrams connected 
and disconnected, with or without insertions of operator O. The logarithm of the partition 
function 

(3.2) 

has the meaning of the generating function for the Green functions of operators O. Dia- 
grammatically, it is given by only connected diagrams with insertions of operator O. In 
the case of many sources, ^[Ji, J 2 , •■•], the corresponding generating function W[Ji, J 2 ,...] 
is given by all connected diagrams with all possible insertions of operators Oi, O 2 , etc. 

Let us consider an operator of the form 


0[A\ = exp 


dx M{x)y[A] , 


(3.3) 


where M{x) is some classical field and T[A] is a composite operator of fields A. According 
to the definition of the quantum average one can consider the vacuum matrix element of 
the operator (3.3) as the partition function evaluated on the “classical sources” M{x), 


{T 0[A]> 


1 


DA M{x)y{x) 


Z[M] 

W 


(3.4) 
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The T-ordering on the left-hand-side of eqn.(3.4) allows us to apply the functional inte¬ 
gration quantization. For brevity, in the following we do not explicitly denote T-ordering, 
assuming it for every quantum average. Exception is made only in sec. 3.3, where we discuss 
the quantum averages of operators without T-ordering. 

The matrix element (3.4) can be presented as an exponent of the generating function 
(3.2). By definition, the generating function can be presented as a series of correlators 

(0[^]) = (3.5) 

= exp^y' dxM{x){y{x)) + ]^ j dxi^ 2 M{xi)M{x 2 ){y{xi)y{x 2 )) 

/ dxi^2,'iyi{xi)M{x2)yi{xy){y{xi)y{x2)y{x^)) + , 

where the factorial coefficients are the symmetry coefficients resulting from the permutation 
symmetry of the correlators. Thus, the exponentiated series for the operator 0[A\ is given 
by only connected diagrams with arbitrary number of operators y convoluted with the 
“classical sources” M{x). 

The fact that VF[M] is given solely by connected diagrams does not imply that the 
original series for (O) contains disconnected diagrams. Both series contain only connected 
diagrams, but connected to different operator vertices. Since, the operators y and O can 
have very different properties, the series of diagrams contributing to (O) can relate to the 
series of diagrams contributing to bF[M] in a very non-trivial way. The case of non-Abelian 
exponentiation, the main subject of this article, is an example of such a non-trivial relation. 

Exponentiation of the perturbative series for an operator of the type (3.3) is the funda¬ 
mental property of the perturbative expansion. It is founded only on the relations between 
the symmetry coefficients of various Eeynman diagrams. In its own turn, the later is the 
consequence of the perturbative approach to functional integration, i.e. the property of 
the expansion of the action exponent around its Gaussian part. It is important to mention 
that the structure of the propagator does not affect the symmetry coefficients and, thus, 
does not inffuence on the exponentiation of the diagrams. 

3.2 Exponentiation of Wilson lines in QED 

Let us demonstrate the application of GF exponentiation of the Wilson line in an Abelian 
gauge theory (QED). 

The starting point of the approach is to present the operator, namely the Wilson line, 
in the form (3.3). In QED the path-ordered exponent (1.2) is equal to the usual exponent, 

= exp dr 7 ^(r)A^( 7 (r))^ . (3.6) 

Gomparing this expression with eqn.(3.3) we find that the operator y{x) is just the field 
A^{x), while the source M(x) is a classical source of radiation, 

M^{x) = ig [ dr j{T) — x). (3.7) 

Jo 
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Evaluating the expression (3.5) with the sources (3.7) we obtain the exponentiated 
expression 


= exp 



dri 



dT 2 Y^ ( ti ) 7^2 {t 2 ){A^^ ( 7 ( 7 - 2 ))> 


(3.8) 



(^Ml(7(7'l))^At2(7(7-2))^/.3(7(7-3))^M4 ( 7 ( 7 - 4 ))) + •••), 


where the correlators with odd number of photons are omitted in consequence of Furry’s 
theorem. Diagrammatically the argument of the exponent is given by connected diagrams 
with an arbitrary number of external photons positioned on the path of the "Wilson line. 

The expression (3.8) is already suitable for the further consideration, however, it can 
be rewritten in the more traditional form. Indeed, using the symmetry of the multi-photon 
correlators under the permutations of fields we combine multiple sources to a single path- 
ordered source, 


= exp - 5 ^^ dny dr27^nn)7''"(7-2)(^/.i(7(n))^/.2(7(7-2))) (3.9) 


1 /■! /•! 4 


+g 



0 JtiJt2Jt2, 


\\dTiY^{Ti)\{A^^{-i{Ti))A^Jcf{T2))A^^{-i{T:i))A^Jci{Ti))) + 


2=1 


The symmetry coefficients of (3.8) are canceled by factorial multipliers resulting from the 
ordering procedure. The expression (3.9) has the common form of the exponent of all com¬ 
pletely connected (connected in the absence of "Wilson line) diagrams with unity symmetry 
coefficients. 

We stress that within GF exponentiation the contour of a Wilson line plays no role. 
Therefore, the exponentiation property holds for any contours, including cusped, self- 
crossed and disconnected. However, for the contours with singularities it is convenient 
to introduce independent sources for every individual smooth segment. This allows one to 
reveal the influence of the contour singularities on the diagrams explicitly. 


3.3 Exponentiation of real exchanges 

Diagrams with real particles exchanges play an important role in applications of quantum 
field theory, especially, in the consideration of final state interactions of hard processes. 
Typically, on the operator level the real exchanges appear as an insertion of the complete 
set of asymptotic states. Alternatively such operators can be given by the (usual) product 
of T-ordered operators 

(T Oi[A] ^ \X){X\T 02[A]) = {{T Oi[A]) (T 02[A])). (3.10) 

On the diagrammatic level it implies that diagrams describing (3.10) contain two parts. 
The interactions inside these parts are presented by the Feynman propagators Xp, while 
the interactions between these two parts are presented by the positive-frequency part of 
the Pauli-Jordan function, 
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We remind that the standard functional integration approach assumes the considera¬ 
tion of only T-ordered operators. The diagrammatic expansion of the functional integral 
can be expressed via the differential reduction exponent (see [12]) 

m 

where Af(x,x') is the Feynman propagator, Sint is the interaction part of the action. In 
the reduction exponent (3.11) the T-ordering of operator is reflected through the Feynman 
propagator. 

In fact, the reduction exponent is the general form of any diagrammatic expansion, 
but not only the T-ordered as it follows from the functional integral. In order to obtain the 
diagrammatic expansion for (3.10) we can adjust every operator to an independent copy of 
a quantum field theory. In this case, the virtual interactions are generated independently 
by copies of the functional integral. The real particle exchanges between these copies can be 
added to diagrammatic by the extra reduction exponent with a real propagator. Formally, 
we have 



0[A\e 


Sir,t[A\ 


A=0 


(3.11) 


((T Oi[^])(r 02[A])) 


Z2[0] 





^1,2=0 


(3.12) 


where we omit the arguments of fields for brevity. In this way, the product of T-ordered 
operators is presented as a T-ordered product of operators in the modihed theory, 


((T Oi[A]) (T 02[A])) = {T Oi[A] 02[^])mod (3.13) 

= [ DA 1 DA 2 0i[^i]02[^2] e^["^i]+^I^2]+/AiA(+)A2^ 

^mod[UJ J 

Such a trick has been applied in ref. [13] for calculation of Drell-Yan soft factor. 

Let us consider two T-ordered operators of the form (3.3). In the course of GF expo¬ 
nentiation we have the chain of equalities 

{{T Oi[A]){T 02[A])) = (T ^ ^w^oAMi,M2\_ ( 3 ^ 4 ) 

The modihed argument of the exponent reads 

ITn,od[Mi,M 2 ] = y'dx(Mi(x)(rTl[^l])mod + M 2 (x)(TT 2 [^ 2 ])mod) (3.15) 

+ J dxi,2(^Mi(a:i)Mi(x2)(rTi[^i]Ti[^i])mod 

+ iM2(xi)M2(x2)(rT2[^2]T2[^2])mod 

+Mi(xi)M 2 (x 2 )(r Tl[^l]T2[^2])mod) + 

where the arguments of operators are omitted for brevity. 
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Any expression in the modified theory can be transformed back to the usual theory 
using the rule (3.13). Therefore, eqn.(3.15) transforms to 

W[Mi,M2] = J dxiMi{x){Tyi) + M2{x){Ty2)) (3.16) 

+ J dxi,2(iMi(xi)Mi(x2)(r TiTi) + ^M2(a:i)M2(x2)(rT23^2) 

+Mi(xi)M2(x2)((r Ti) {T T 2 ))) + 

where the arguments of operators are omitted for brevity. We conclude that the diagrams 
with real exchanges exponentiate in the same way as the usual virtual diagrams. 

4 Exponentiation of Wilson lines in non-Abelian gauge theories 

The exponentiation in non-Abelian gauge theories is a more involved problem than the 
exponentiation in Abelian ones. In comparison to QED there are two sources of compli¬ 
cation. The first complication comes from the involved form of operators y. The second 
complication comes from the (color-)matrix structure of the non-Abelian Wilson line. The 
structure of operators y and the generating function was elaborated in ref. [9]. However, 
in ref. [9] the matrix issues of the non-Abelian exponentiation were missed. In this section 
we present a detailed derivation of the non-Abelian exponentiation. 

We remind that within the framework of GF exponentiation the only important in¬ 
formation is the path-order of gauge fields, but not their coordinates. The term Wilson 
line, which is used all over the article, denotes rather a Wilson curve on arbitrary path, 
than a straight Wilson line, as it may be suggested from the term. Generally, the path can 
be arbitrary difficult, or even non-analytical. However, we suppose that every individual 
Wilson line is smooth. 


4.1 Wilson line as exponent 

The starting point of GF exponentiation is to present the Wilson line (1.2) in the form 
(3.3). For this purpose we use the following exponential representation for the Wilson line. 



^ ''If] 

where A = A^{'^{Ti)) dxi and tq = r, with A^{x) = fa A“(x). The operator ad^i is 

dehned as ad^iA = [A, X]. The detailed derivation of this relation can be found in ref. [14] 
The representation (4.1) reveals several important properties. The main of them is 
that the operator y consists only of completely nested commutators of generators of gauge 
fields. It leads to the color-connectivity of webs, the property that has been the defining 
property of the webs for many years starting from [3]. 
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The completely nested commutator structure of the operators also implies universality 

of the representation (4.1) for Wilson lines of different group representation. Indeed, the 

f [ fl \ f] 

commutator of generators is proportional to the structure constant [ta ,^ ] = ifabcic ■ The 
structure constant is independent of group representation. Therefore, the representation 
dependence of the operator y is concentrated in a single generator, which can be moved 
out of the operator. The equation (4.1) takes the form 


= exp 


(4.2) 


where are operators independent on the representation of the Wilson line. Here the 
definition of the operator V is slightly different from the corresponding definition in ref. [9] , 
where the corresponding operator was defined with the integration over tq removed. The 
representation (4.2) also has an advantage that all the matrix structure of the Wilson line 

jfi 

is concentrated in the single generator ta ■ 

In the case of multiple Wilson lines, it is convenient to merge all Wilson lines to a 
single exponent. With this purpose we demand that every separate Wilson line acts in 
a separate matrix space. Then the whole composition of Wilson lines is a matrix of the 
reducible representation /i (8i /2 <8) ... <8> /at, where fi is the representation of z’th Wilson 
line. This trick is often used to simplify the consideration of web diagrams, see e.g.[5, 7]. 
In order to simplify the notation we introduce the generator in this space with appropriate 
numbering 


Ta = 


Ah] 

1 (8) t 


[h] 


1, A = 1, ...,dimG, 

1 , A = dime + 1,..., 2dime, 


(4.3) 


1 (g) 1 (g) .. (g) A = {N — l)dimG + 1,..., A^dime 


where dime is the dimension of the gauge group or, equivalently, the number of generators. 
We adopt the convention that the labels in the joined space are capitalized, while the labels 
in the irreducible spaces are denoted by the lowercase letters. Simultaneously, the path 
dependence of the operator I4y“ is adjusted to the corresponding sector of label A, 




V^, A=l, ...,dimG, 

V^, A = dime + 1,..., 2dimG, 


H = (iV - l)dimG + 1,..., A^dim^ . 


(4.4) 


In this representation the product of several Wilson lines reads 

. xM x...x(^%)^ =^jj={expTAV^)^j, (4.5) 

\ / liji \ / 12J2 ^ ' iNjN 

where I and J are indices of the joined matrix space. 

The representation (4.5) resembles the required form (3.3), where generators Ta play 
the role of the matrix sources for operators V^. For the general discussion of GF exponen¬ 
tiation it is convenient to keep the infinite sum of operators (4.1) as a single object V^. 
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While for the perturbative analysis and practical applications it is convenient to split the 
operator on individual terms of a fixed perturbative order, 

CX) 

= ( 4 . 6 ) 

n=l 

where Vn ^ g^- The expressions for the first few operators Vn, as well as, Feynman rules 
for them are presented in the appendix A. 

4.2 Matrix exponentiation 

According to the general discussion of sec. 3.1, the connected diagrams with all possible in¬ 
sertions of the operators are the only diagrams contributing to the exponent. However, 
in the case of non-Abelian gauge theory it is not entirely correct. The point is that in the 
non-Abelian gauge theories the operators in diagrams are caused by the matrix sources T^. 
This complication spoils the usual relations between connected and disconnected diagrams, 
and prevent the straightforward exponentiation. 

The direct way to bring the diagrammatic series into the exponentiated form is to 
consider its logarithm. The perturbative expansion of the logarithm (of the perturbative 
series) mixes up the diagrams such that it is very difficult to find out on a general level, 
which parts of diagrams cancel and which do not. Therefore, this way necessarily leads us 
to the analysis of individual Feynman diagrams, the procedure which that we try to avoid. 
This approach was elaborated in ref. [7], and indeed it appears to be not efficient. 

In order to solve the matrix complications of the non-Abelian exponentiation in the 
most efficient way, we split the consideration onto two independent tasks. First, we consider 
the exponentiation of a scalar version of the non-Abelian Wilson line. That grants us the 
most promising starting point for the exponentiated diagrammatic. Second, we generalize 
the scalar version of Wilson line and its exponentiated expression on the matrix form. In 
this way we obtain the complete matrix exponentiated expression without lost of efficiency 
of GF exponentiation. 

The first point of our program is to consider a scalar operator 

(j) = exp {M^Va) , (4.7) 

where is a scalar source. We emphasize that the term scalar in this section denotes 
an object, which does not carry matrix indices (/J), although the object can be a vector, 
or a tensor in the color space. For example, the operator in eqn.(4.2) is a scalar, while 
the generator is a matrix. 

The operator cj) has many common properties with the non-Abelian Wilson line <h. So, 
one may say that is a scalar image of the non-Abelian Wilson line. According to the 
discussion of sec. 3.1 the average of (f) can be presented in the form of exponent 

(</-) = §|^ = expW[M]. (4.8) 

The function VF[M] is given by all connected diagrams with insertions of operator Va- The 
exponentiation of the scalar image of Wilson lines is, as simple as, the exponentiation of 
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Wilson lines in the Abelian gauge theory, discussed in sec. 3.2. The only difference is that in 
the non-Abelian gauge theory the operator V is given by infinite series of operators (4.6). 

In order to fulfill the second point of our program, we need the formal definition of the 
matrix generalization procedure for a function of several arguments. It can be presented 
as action of the matrix shift operator on the scalar image. 


, rpA d 

f{T)ij = {e^ ^ 


1-3 


x=0 


(4.9) 


It is straightforward to check that right-hand-side of eqn.(4.9) satisfies all standard demands 
on the matrix function f{T). 

The Wilson line ‘h is obtained from its scalar image (4.7) with the help of operation 
(4.9) 


/ 

\ j. 

1 e 

4> 

J IJ 


= 

Therefore, the average of Wilson line is given by 

Z[M] 


M=0 


^Ij) = 


rpA 5 

e Tm^ 


IJ Z[0] 


M=0 


Zij[T] 

Z[0] ^ 


(4.10) 


(4.11) 


where Z[T] is the matrix generalization of the partition function by means of the procedure 
(4.9). 

As one can see from eqn.(4.9), the generalization of the partition function to ma¬ 
trix sources is rather straightforward procedure. However, the relations between diagrams 
within Z[T] are not the same as within Z[M]. This happens due to the symmetrization of 
matrix variables by the shift operator (4.9). As a result, the originally disconnected dia¬ 
grams are entangled by their matrix structure (see explicit example in eqn.(4.24)). There¬ 
fore, the logarithm of the matrix partition function Z[T] is not the matrix generalization 
of the function IT[M] (which we denote as IT[T]). In symbolic notations this statement 
reads 


Zij[T] 

Z[0] 


1 


pA S 


I e 

Z[0] V ) IJ 


W[M\ 


M=0 




.W[T]\ 


' IJ 


where 


Wij[T] = e 


IJ 


W[M] 


M=0 


(4.12) 


(4.13) 


The inequality (4.12) has been overlooked in ref. [9], therefore, the final result of ref. [9] 
presented there is incomplete. 

Considering the left-hand-side of eqn.(4.12) we conclude that although the function W 
is not the complete result of the exponentiation, but still it is the only function that can 
appear in the argument of the exponent. Taking into account that the leading term of the 
perturbative expansion is W, we can present the left-hand-side of eqn.(4.12) as 


ZijjT] 

( f,W[T]+5WlT]\ 

Z[0] \ Jij- 


(4.14) 


- 12 - 












where 5W is a function of W. 

One can see that the function W plays an exceptional role in the GF exponentiation 
procedure. For this reason we call W the kernel of matrix exponentiation (MEK). 

The correction terms 6W we name the defect of matrix exponentiation (for shortness 
we often call it the defect). Mathematically, the defect arises from the reordering of color- 
matrices during the formation of matrix exponent. Accordingly, the structure of the defect 
resembles the famous tail of commutators in Baker-Campbell-Hausdorff (BCH) formula. 
The defect is an algebraic function of VF[r], in the same manner as the higher terms of 
BCH series are algebraic functions of the previous ones. The formal definition of the defect 
reads 


5Wij[T] 



M=0 




Ine^I^] 

IJ 


M=0 


(4.15) 


Thus, we have shown that the (T-)product of Wilson lines can be presented as a matrix 
exponent 

( H). . - (-i-ft). , > = (-1-/4 = . ( 4 ^ 16 ) 

^ ' HJl ^ ' i-nJn ^ ' IJ 


where functions W and 6W are defined in (4.13) and (4.15) respectively. This expression 
is the main result of the article. In the following we discuss the properties of MEK and the 
defect, and demonstrate the relation between the presented approach and other approaches. 

The form of the exponent (4.16) would necessarily appear during the exponentiation 
of any matrix object. The most important point of exponentiation is the selection of 
an efficient leading term, which is MEK in GF exponentiation. As we already told, in 
ref. [7] the whole perturbative series has been chosen as a leading term, and this choice 
leads to unnecessary complications. Our choice of MEK is inspired by many remarkable 
features of W. Above all, MEK is the generating function for webs in the case of the 
scalar image of Wilson line. It implies that MEK already contains essential features of 
the exponentiated expression, except the matrix issues. Another example of important 
features favoring MEK is that MEK contains all color connected diagrams (see details in 
ref. [9]). The color connectivity is a defining attribute of web diagrams [8]. All in all, MEK 
is the best candidate for the leading term and, as we demonstrate later, such a conjecture 
is well-founded. 

The structure presented in eqn. (4.16) is can be also traced in the replica exponen¬ 
tiation [5, 6, 8, 10]. Within the replica exponentiation the diagrams contributing to the 
exponent always has two contributions: the part proportional to unity in the replica space, 
and the rest, which has more involved structure (for explicit examples see e.g. equations 
(50),(61),(63) and (65) in [8]). The part proportional unity, after the summation over 
replica indices, is linear in Arep, and therefore, directly contributes to the exponent. In 
GF exponentiation this part is given by MEK. The rest is an arbitrary polynomial of Arep, 
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the linear part of which contributes to the exponent. In GF exponentiation this contribu¬ 
tion is given by the defect. We emphasize that in contrast to the replica exponentiation, 
in GF exponentiation the defect (4.15) is given by explicit expression, that significantly 
simplify the consideration of webs. The detailed comparison of the approaches, as well as, 
an example-calculation are given in sec. 5. 


4.3 Structure of MEK 

MEK is a matrix generalization of the function VF[M] (4.8). Using the explicit expression 
for generating function (3.5) we obtain 

(rTiArriB\ (rjiArr'BrpC\ 

= Tfj{V^) + 4- —h^iyAyB'^ ^ i^- hl(^V^yByC^ ^ ^ 

^. o I 

where we have used that the correlators are symmetric over permutations of operators. 

Important to note that in the case of real gluon exchanges the matrix generalization of 
generating function is more involved. The product of T-ordered operators is not symmetric 
under permutations. Therefore, the expression (4.17) would contain terms with explicit 
symmetrization of generators. 

Let us pass from the joined notation for multiple Wilson lines (4.3-4.4) to the consid¬ 
eration of individual contributions. We continue to assume that generators related to a 

r f 1 

separate Wilson line act in the separate matrix spaces. So, the generator ta is a matrix 
acting in the space {ik,jk) and, MEK is a multimatrix in spaces {ik,jk) for A: = 1 to A^. For 
brevity we continue to denote the multi-index by capital letters, i.e. = Wjj. 

We have 

N 

= (4.18) 

7 ' '^kjk 

k=l 



where the missed indices on the right-hand-side of the expression should be given by the 
unity matrices. Also we will omit the representation indication [/] on generators, since it 
is hxed by denoting the matrix indices. 

In the case of multiple Wilson lines, MEK reads 


N 


'rH = E‘L.(U“)+E 


N ,a jh 
^hjk 


k=l /c,/=l 

k^l 

N .a j-b J.C 

I ^ ^ _jkjk '^ijl '^mjm 


2 ! 


N Ua^b\ 

ya yfe \ I J-llML / yo- yb 


fc =0 


k,l,m=l 
k^l^m 

k=l 


3! 


N fC 

kKv‘\ + E r "E utUu;) U 19 ) 


k,l=l 

k^l 


2 ! 


iMh / yb yc\ , 
31 \ Ik 'ik 'ik / ^ 
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Figure 1. Diagrams contributing to MEK for three Wilson lines at 0{g^) order that connect all 
three Wilson lines (diagrams with permutations of Wilson lines should be added, as well as, web 
diagrams of 0{g'^) order with a virtual loop). Blobs with Vn denote the vertices (A.l), while the 
empty blob denotes all possible four-gluon tree interactions. Ovals with a “plus” sign denote the 
symmetrization of the vertices. 


This can be simplified by the fixation of the order of individual Wilson line within correlator, 
N N N 


‘Uj = + E + e*S(k,:<) 

k=l k,l=l k=l 

k<l 

N N 


k,l,m=l 

k<l<m 

N 


k,l=l 

k<l 


N 


^ \ ''ll '^'11 / ^ ''ikjk \ Ik '^ik '^ik ^ 


k,l=l 

k<l 


k=l 


where denotes the product of generators symmetrized over the indices 


^{ai...a„} _ J_ 

n! ^ 




cr=perm[ai...an] 


(4.20) 


(4.21) 


According to eqn.(4.20) one has a simple diagrammatic rule: MEK consists of all 
diagrams connecting arbitrary number of operators V located on Wilson lines with unity 
symmetry coefficients. Every operator should be convoluted with the gauge group generator 
of the corresponded Wilson line. If there are several operators adjusted to the same Wilson 
line they should be convoluted with the symmetrized product of generators. 

Eor the actual calculation one inserts the perturbative expansion for the operator V 
in terms of operators (4.6). The operator Vn radiates exactly n gluons. Nonetheless, the 
general diagrammatic rules are still the same: the gluons could be contracted to other ver¬ 
tices, interact with each other, or be contracted to the same vertex, but all vertices V must 
be connected together. An example of diagrams with three Wilson lines is demonstrated 
in fig.l. 
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4.4 Structure of the defect 


The formal definition of the defect is given in eqn. (4.15). However, the form (4.15) is not 
transparent and is hardy applicable. In this section we reorganize eqn. (4.15) in simpler 
form, suitable for calculations. 

In order to simplify eqn. (4.15) we expand the exponent of the generating function in 
a series 


^W[M] ^ I I m + ^ H/g + SITaH^i + ITf ^ 


3! 


OO 

E E 

n=0 {fc} 

J2ii-ki=n 






wh 


(l!)^iA:i! (2!)^2 /j 2! (re!)^"A:,i!’ 


(4.22) 


where we have introduced the short notation 




(4.23) 


The action of the matrix shift exponent (4.9) replaces all sources by the generators 
in the completely symmetric way. This operation mixes the matrices, such that the result 
is not the matrix product of Wn’s. For example, the first term with mixture of matrices in 
eqn.(4.22) is 


/ rpA _^ 

(e 5m- 


IT 1 IT 2 

IJ 


M=0 





^rpBrj^C _j_ rj^ArpC rpB _|_ rpC rj^ArpB'^ A'^'B 


(4.24) 


where the symmetry of correlators is taken into account. One can see that this expression 
is not equivalent to (VT 1 IT 2 + IF 2 VTi)/ 2 , what is needed for the expansion of exp(iy), but 
contains extra terms. The extra terms are proportional to commutators of generators and 
give rise to the defect of exponentiation. 

Let us introduce the special notation for symmetrized products of the form (4.24). We 
denote 


{Wn]lJ 

{WkWn-k}lJ 

{WkWiWn-k-l]lJ 


^IJ 

■ • *-^71} 

^IJ 

^IJ 


( .. 1/^'=) (I/^'=+L.. ), 

( .. 14^'=) ( .. 1/^'=+* ){V^ 


fc + Z + 1 



etc. 


(4.25) 


where denotes the symmetric product of generators, see (4.21). 

The operation (4.25) does not have a special name, nonetheless it often appears in 
course of diagrammatic resummations and related problems. In ref. [7] this operation serves 
as a generalized product of diagrams needed for the diagrammatic exponentiation. Authors 
of ref. [7] shows that the result of the symmetrized product (4.25) has the meaning of the 
sum of all joined diagrams with all mutual ordering of gluons attached to Wilson lines. 
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In the notation (4.25), the matrix generalization of the partition function reads 


Zij[T] = (e 


pA 5 


,W[M] 


ij M=o (l!)^i (2!)*^2...(n!)^’*/ci!A:2!...A:„! 

n=0 {A:} 

Y,ii-ki=n 


E E 


. (4.26) 


In the same notation MEK and the partition function (4.26) are 

{W[Tf}ij 


iy„|Tl = E^. 

k=l 


Zij[T] = Y.- 


k=0 


k\ 


(4.27) 


where the operation {..} for {IE[T]^} should be taken individually for every term of per¬ 
turbative expansion. 

Finally, we evaluate (4.15) and obtain the expression for the defect. For the practical 
applications it is convenient to present the defect as the series 


5Wij[T] = Y,^kWij[Tl (4.28) 

k=2 


where 6kW ^ W^. 


The first few terms are 


62 Wij[T] 

w/jm 

SaWjAT] 


, {{wnw)_^^{w{wA 

6 ^ 4 


IJ 


3 V /IJ 


(4.29) 


m.j ({tr»}tr)„ + (ir{irn)„ 


4! 


12 


+ 


w'AW^) + {w^{wA) + {w{wAw^ (w- 


ij 


The equations (4.29) can be significantly simplified by means of recursion. With this 
aim we consider eqns.(4.29) as equations on the totally symmetrized products {W^}, as 
functions of SW. Substituting the solution again into eqns.(4.29) we found a simple recur¬ 
sive relation for the defect 


1 1 

&ulTI = - E t! E (Pr-Cir)^^ , (4.30) 

k=2 i>l 
E*=»4 

where 5iW = W. In the next sections we present explicit calculation of the defect for 
particular diagrams. 

In fact, the matrix structure of the defect 5nW is proportional to (at least) (n — 1) 
commutators of generators. It follows from the recursive equation (4.30). Therefore, if 
matrices T commute, all term of the defect are zero. In this way many contributions to 
the defect are zero, especially in for diagrams that involve many Wilson lines. 
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Due to the property that the defect is proportional to the commutators of generators, 
the color factors appearing in the defect are color connected. Together with the statement 
that MEK is color connected, see discussion in [9], it proves the non-Abelian exponentiation 
theorem formulated in [8]. 

It is interesting to consider the gauge invariance of the exponentiated expression. Gen¬ 
erally, the correlator (4.16) is not gauge invariant, and therefore, there are no special prop¬ 
erties for MEK or for the defect. However, if one considers a color singlet configuration, 
or the open color indices located at spatial infinities, the correlator (4.16) is perturba¬ 
tive gauge invariant (for the later case, only in non-singular gauges, see e.g. discussion in 
ref. [15]). Consequently, the exponentiated expression for such a gauge-invariant configu¬ 
ration, is also gauge invariant. Indeed, at order a the statement is obvious. Given that, 
the gauge invariance for higher orders can be proven by iterations. In its own turn, the 
gauge invariance of the exponentiated expression leads to gauge invariance of the defect 
and MEK independently. 

Let us sketch the proof of independent gauge invariance for MEK and the defect (given 
that the correlator under consideration is perturbative gauge invariant). The proof is made 
by interactions, comparing order-by-order the (gauge invariant) exponentiated expression 
and 1K[T] -|- (51K[r]. Comparing the leading terms we observe that the leading term of 
MEK is gauge invariant. The next-to-leading term is given by the sum of next-to-leading 
term of MEK and leading term of the defect. The leading term of the defect is gauge 
invariant, because it is composed from the leading terms of MEK, see (4.30). Therefore, 
the next-to-leading term of MEK is gauge invariant. Then one can construct iteration and 
prove that MEK and the defect are independently gauge invariant at any order. 

5 Application of GF exponentiation 

In this section we compare the GE exponentiation with the diagrammatic and the replica 
exponentiations. The examples presented in the section clarify the role of the defect and 
MEK in the exponentiation procedure. 

In order to compare approaches we consider two half-infinite Wilson lines, i.e. the cusp. 
The cusp is a popular playground for the testing of various methods. Recently, it has been 
evaluated up to the three loop order in the general kinematics in QCD [16]. The evaluation 
of exponentiated diagrams at one and two-loop order can be found in many articles, e.g. in 
refs. [17-20] one can find all necessary details and the explicit expressions. In the course 
of comparison we do not evaluate loop-integrals. Instead, we compare the unintegrated 
expressions for the diagrams and show the relations between different approaches. 

In order to demonstrate the effectiveness of GE exponentiation we consider the spe¬ 
cial class of diagrams called the multiple gluon exchange webs (MGEW). These are the 
diagrams with neglected interaction between gluons. Such an approximation is simplest 
but non-trivial part of the exponentiated expression. Particulary, the simple structure of 
MGEWs results to the possibility to evaluate MGEW-loop-integrals at high orders without 
significant efforts, see e.g. [21, 22]. MGEWs were studied in details within the replica expo¬ 
nentiation in refs. [11, 22]. For simplicity, we consider the configuration of light-like Wilson 
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Figure 2. The web diagrams contributing to the correlator of two Wilson lines. The diagram B 
should be taken with the modified color factor. 



Figure 3. The diagrams contributing to MEK for the cusp. The diagram B does not equal the 
diagram B in fig. 2. Diagrams C, D and E, as well as, the diagram A are equal to the corresponding 
diagrams in fig. 2. Diagrams F, G, H and I are zero due to the convolution of antisymmetric three- 
gluon vertex with the symmetrized product of generators. 


lines, and show that in such kinematics MEK is given solely by one-loop diagrams. There¬ 
fore, all the higher order contribution are generated by the defect, and can be evaluated 
by algebraic manipulations. 


5.1 Cusp at two loops 

Let us consider a cusped Wilson line, which approaches from infinity to the origin along 
vector ui, and then continues to infinity along vector V 2 - For definiteness, we consider the 
Wilson line in the fundamental representation of SU{Nc). We denote 


^ .. ) = SijT{vi,V 2 ) = dij exp [CfE{vi,V 2 )] , (5.1) 


where is the half-infinite Wilson line pointed from the origin along the vector v. For 
simplicity, we neglect self-interaction of Wilson lines, and use the Feynman gauge. 

The webs contributing to E within the diagrammatic exponentiation [3, 4] are pre¬ 
sented on fig. 2. The webs are to be equipped with modified color factors (2.2). At this 
order of accuracy, the only diagram, whose color factor differs from the original, is diagram 
B. Its color factor is —CfCaI‘^- 

In order to apply GF exponentiation as it is described in sec. 3, we consider both Wilson 
lines in separate matrix spaces, see sec.4.1. At two-loop accuracy, only one combination of 
generators appears in the exponent. We denote 
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where dots stay for other compositions of color indices arise at three-loop order (see example 
in the next section). Here, the exponent is the matrix exponent, where generators tj^j^and 
ti^j^ act in different matrix spaces. 

The label C (as if charge-conjugation) on the first generator in eqn.(5.2) recalls that 
these generators should be multiplied in the reverse order, from right to left. Such a 
prescription is needed to match out definition of the operator Hy, which utilizes the path¬ 
ordering from the origin, but not from inhnity. In other words, the proper color and path 
ordering for conjugated Wilson line is = exp (—where the minus sign is the result 
of the conjugation of operator V. 

Contracting eqn.(5.2) with and using that function is proportional to 6°’^, we 
obtain the relation (at this order of accuracy) 

E<^^(vi,V2) = S^^E(vi,V2). (5.3) 

In its own turn, the function E°‘^ is given by the sum of MEK and the defect (4.16). 

Let us remark, that one can consider the cusped Wilson line in another way. Instead 
of consideration of every Wilson line in separate color space, as it is done in (5.2), one can 
consider a single Wilson line, whose path-dependent part is split, = 14“ -|- 1^“. Within 
such consideration, one can keep a single generator in (4.2) by the price of two operators 
V. These two approaches are equivalent to each other, although diagrams would distribute 
differently within the perturbative expansion of MEK (4.20). 

In the following we perform the diagram-by-diagram comparison for both sides of (5.3). 
The diagrams contributing to MEK are presented in fig. 3. In the context of eqn.(4.20), the 
diagrams from A to E are produced by the perturbative expansion of the double correlator, 
while the diagrams F and G are produced by the triple correlator. Contribution of the 
average of an operator V is zero in the absence of self-interaction. 

Let us compare diagrams one-by-one. The comparison is simpler to perform in the 
position representation. The gluon propagator reads 

Kb y) = y) = ^ _ y)2 + ^ 0 ) 1 -^ ’ 

where e is the parameter of the dimension regularization, d = 4 — 2e. The expressions for 
the vertices Vn are presented in appendix A. 

The diagrams A and E stay the same in both approaches. For the following consider¬ 
ation we need the one-loop expression for MEK, 

roo 

Wi.ioop = K)fijKi 2 j 2 ^iKp =/ dxdy A{vix,V 2 y), (5.5) 

Jo 

where Vij = {vi ■ Vj). 

The diagrams C and D are also equal in both approaches. Let us demonstrate it 
explicitly. The diagram C in fig. 2 reads 

d^jOpEc = {igf (5.6) 

\ / ij I 
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The kinematical part equals to 

yc = 


poo r 

J dydxi ^2 J d^^z V^up{z)v^v'(v^A{z,viXi)A{z,viX 2 )A{z,V 2 y) 6 {x 2 < xi), 

where V^^^{z) is the Feynman rule for the triple-gluon vertex. The diagram C in fig. 3 
reads 


~ab _ 1 (* 5 -) (* 5 ^ 


Tpao _ 

— 


^ jadc jacd^j ^ jcdbj^^ ^ab 




(5.7) 


where we have used the Feynman rules for V 2 (A.4) and the symmetry properties of the 
three-gluon vertex. The factor 1/2 in front of expression (5.7) is the symmetry coefficient 
of the diagram. Therefore, we find that = d°‘^Ec- The consideration of the diagram D 
is analogous. 

The diagrams F, G, H and I have not analogues in the diagrammatic exponentiation. 
These diagrams equal zero, due to the convolution of the symmetric combination of gen¬ 
erators with the anti-symmetric color structure of the three-gluon vertex, or V 2 vertex. 
Explicitly, we have 

= 0, 


*Ui 




(5.8) 


where the kinematical part is 

poo p 

-Af = y dydxi^2 J d^^z Vpi,p{z)v'^ViV^A{z,viXi)A{z,viX2)A{z,V2y). 

The consideration of diagrams G, H and I is analogous. 

Finally, the diagrams B in fig. 2 and in fig. 3 are not equal. This is the manifestation of 
the matrix origin of Wilson lines. Within the diagrammatic exponentiation the diagram B 
reads 


Eb = 


(5.9) 


poo 

X / dxi^2dyi,2^{viXi,V2yi)A{viX2,V2y2)0{xi < X2)0{y2 < yi). 
Jo 

While, in GF exponentiation the contribution of the diagram B to MEK reads 


■pab 

^S,MEK 


-1 


• 2 \ 2 


= ^(^1 


poo 

/ dxi^2dyi,2 
Jo 


(5.10) 


xA{xi,yi)A{x 2 ,y 2 ) i0{xi < X 2 ) - 9{x2 < xi)) {9{yi < ^ 2 ) - 9{y2 < yi)) ■ 

The factor 1/2 in front of the expression is the symmetry coefficient. 

Expressions (5.9) and (5.10) differ significantly. The difference is the defect of expo¬ 
nentiation, which comes from the squaring of the one-loop contribution (5.5). Evaluating 
the expression (4.29) we obtain 


<52^ = ^! 


{^gr 




{eE + {eE + EE)i. 


1232 




aj.6\ 


)^232 




(5.11) 


poo 

/ dxi^2dyi,2A{viXi,V2yi)A{viX2, ^ 2 ^ 2 )- 
Jo 
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Simplifying the color structure we find 

^l^defect = (ixi,2dyi,2A(uiXi, U2yi) A(uiX2, ^27/2)■ (5.12) 

Adding the defect contribution (5.12) to MEK contribution (5.10) and simplifying the 
integrands we obtain 

+ -^S^defect = (5.13) 

Therefore, the sum of diagram B for MEK and the defect reproduces the diagram B within 
the diagrammatic exponentiation. 

Thus, we observe the complete agreement between GE exponentiation and the dia¬ 
grammatic exponentiation. We stress that the correspondence between approaches is not 
trivial. Namely, the diagrams contributing to MEK have an unique kinematical part. Only 
in the combination with the defect, the standard diagrammatic (2.1) is restored. It implies 
that GE exponentiation exponentiates not only the color part but also the kinematic part. 
So, the method separates unique contributions of two-loop integrals, from the powers of 
one-loop contribution. 

It is also instructive to compare GE exponentiation with the replica exponentiation. 
The diagrams appearing in replicated theory are the same as shown in fig. 3. The expressions 
to compare in the form of eqn.(5.2) can be found in ref.[8]. 

The contributions of diagrams A, C, D, E, F and G equal in both approaches. The 
equivalence happens because these diagrams are irreducible in the absence of Wilson lines. 
Therefore, only the fields of the same replica quantum number are presented in the diagram. 
In other words, these diagrams are proportional to unity in the replica space, that after 
evaluation of the diagrams results in the common factor Wep (number of replicas). Effective 
vertices at A^rep = 1 are the same in both approaches (compare our definition of V in (4.2) 
with the dehnition (20) in ref. [8]. Also compare explicit expressions for the vertices Vn 
(A.l) with (22a-22c) in ref. [8]). Thus, the rest parts of diagrams are identically equal in 
both approaches. 

Diagrams H and I are zero in both approaches due to the convolution of symmetric 
and anti-symmetric color factors. In ref. [8] this fact is demonstrated in eqn.(48). 

Diagram B in the replica exponentiation is presented in eqn.(52) of ref. [8]. Using our 
notation this equation reads 

(',-„ 2\2 roo 

EB,rep = - j dxi^2dyi,2^{viXi,V2yi)/^{viX2,V2y2) (5.14) 

X [Wep {0{xi < X 2 ) - 9{x2 < Xi)) {9{yi < 2 / 2 ) - 9{y2 < yi)) + Nj-ep (N^ep - 1)], 

where we have taken into account that in ref. [8] both Wilson lines are incoming. Gomparing 
this expression with (5.10) and (5.12) we observe that 

E'B,rep — EtrepE’^ypBK ~ Arep(Arep “ l).E'B^defect■ (5.15) 

Therefore, the contribution of MEK and the defect are clearly distinguished within the 
replica exponentiation. Namely, the diagrams with only a single copy of field (i.e. with 
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Figure 4. MGEWs contributing to the cusp at three-loop level order. The symmetric diagrams 
are not shown. For light-like Wilson lines only the first diagram is not zero. 

all internal vertices proportional to 6nm in replica space) are the contribution to MEK. 
These diagrams are proportional to the first power of A^rep- The rest contributions are 
proportional to a polynomial in Wrep and represent the contribution of the defect. The 
polynomial in A^rec reflects the reordering of color matrices during exponentiation, which 
is taken into account explicitly by definition (4.15) in the case of GF exponentiation. 

5.2 Multiple gluon exchange ^vebs for light-like Wilson lines 

The webs, as well as, diagrams contributing to MEK are naturally split on subclasses by 
the number of connected parts of the internal graph. The division by this principle is not 
gauge invariant, but it is very effective from the practical point of view. 

The limiting case of such a division is the completely connected diagrams, i.e. diagrams 
that are connected in the absence of both Wilson lines and vertices Vn (e.g. these are 
diagrams A,C,D,E,F,G in fig. 3, and diagram A in fig.l). These diagrams are not influenced 
by the matrix structure of the non-Abelian gauge theory, and do not mix with the defect. 
For this subclass, the diagrams of MEK are in one-to-one correspondence with the diagrams 
of the diagrammatic exponentiation. The limit of completely connected graphs resembles 
the Abelian exponentiation procedure. 

The opposite limiting case of the division is the diagrams without any internal inter¬ 
action of gluons (e.g. these are diagrams B,C,E, and F in fig.l). These diagrams have been 
studied in details in the recent publications [11, 22], and are called multiple gluon exchange 
webs (MGEW). Within GF exponentiation MGEWs are given by the sum of MEK and the 
defect contributions. Moreover, the contribution of the defect is the strongest in compari¬ 
son to the more connected diagrams. In some sense, MGEWs contain the smallest amount 
of a new (in comparison to the previous order of perturbative calculation) information. 

MGEWs are the best demonstration of the effectiveness of GF exponentiation. In the 
case of light-like Wilson lines MEK can be evaluated exactly, and directly contributes only 
at one-loop order. Therefore, for light-like Wilson lines the higher-loop orders of MGEWs 
are given entirely by the defect of exponentiation. 

To start with, we consider the cusp configuration (5.2) at vf = V 2 = 0. MGEWs 
contributing to the cusp up to three-loop level order are shown in fig. 4. The one-loop 
contribution to MEK is given by (5.5). For the light-like kinematic we have 



(5.16) 
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where exponent performs the infrared regularization. For future convenience we rewrite 
the one-loop expression in the form 


Ef = 


Ki 


(27r)i-^ V-^V 


(5.17) 


The two-loop contribution to MEK is given in eqn.(5.10). In the light-like kinematics 
it reads 


jpab 


2 


, 2 \ 2 


j;acd fbcd„.2e^ 

J J ^12 16vr4-2. 


POO 

dxi^2dyi,2 


(5.18) 


\ — / - ^ u 

{9{xi < X 2 ) - 9{x2 < xi)) {e{yi < ^ 2 ) 




The zO prescription does not influence the integral, because the ultraviolet singularity at 
xy —>• 0 is regulated by the dimension regularization. We conclude that the propagator 
attached to light-like Wilson lines can be split on individual components, A ~ 

Therefore, the loop integral is zero due to the antisymmetry of the integrand under the 
permutation of xi and X 2 , or yi and y 2 - 

Thus, in the light-like kinematics of Wilson line at two-loop MGEW contribution is 
given solely by the defect contribution. Evaluating integral (5.12) we obtain 

-^2Uoop,MGEW = • (5.19) 

This result coincides with the calculation performed in the diagrammatic exponentiation, 
see e.g.[20]. 

In the similar way one can show that the contribution of the three- and higher-loop 
diagrams to MEK are zero. Indeed, all the propagators for MGEW with light-like Wilson 
lines are just the product of components, A ~ Therefore, loop-integrals can be 

split on two parts: one depending on the variables related to Wilson line and another 
depending on the variables related to ^V 2 - The interaction with Wilson lines are given by 
vertices Vn- Schematically, the contribution of some n-loop diagram reads 


/■c 

y (color factor) / 




Vk{xi,..)...Vi{..,Xn) 




Xi 


■ '^n 
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dyi,...,, 


Tr(yi) ■■)---Vs{--,yn) 




■ yi " 


(5.20) 


where the sets of vertices {k,..,l) and (r, ..,s) are defined by the diagram, k + ... + I = 
r + ... + I = n, and the summation runs over all independent color structures. According 
to the diagrammatic rules for MEK, there are necessarily several vertices Vn with n > 1. 
The kinematical part of vertices Vn{xi, ..,Xn) is antisymmetric over the permutation of 
arguments. Therefore, the integrals on the right-hand-side of (5.20) are zero. The only 
exception is the case when all vertices are Vi, which is the one-loop contribution (5.17). 

We conclude that in the MGEW approximation, MEK is given solely by the one-loop 
diagram, and reads 


Wlj {vi,V2) 




MGEW 


(5.21) 
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Afc(10-33Ar^^) 


4608 

768 

4608 

5 

-35+93N^-57N^ 

{l-N^){23N'^-22) 

Af2(71-114Af2) 

28800 

23040 

57600 


Table 1. The coefficients for expansion of the light-like cusp amplitude in MGEW approximation 
(5.22-5.23) for SU(A^c) gauge theory. 


where Ki is given in (5.17). 

The fact that at MGEW approximation MEK is given by solely by one-loop diagram, 
leads to the surprising conclusion. Namely, MGEW approximation is gauge invariant (in 
non-singular gauges) in the case of light-like kinematics, as long as, Wilson lines are color 
connected at the origin. Indeed, MEK is given by a single gauge invariant diagram. As a 
consequence, the defect is also gauge invariant, see discussion in the end of sec. 4.4. The 
gauge invariance of MGEW approximation holds only for light-like Wilson lines. Eor the 
Wilson lines off-light-cone the MEK has contributions of higher orders, which are composed 
from different diagrams, and therefore, MGEW diagrams are not gauge invariant. 

The contribution of the defect can be obtained by straightforward evaluation of ex¬ 
pressions (4.29) with MEK (5.21). Eor an arbitrary gauge group the evaluation involves 
the elaboration of higher Gasimir operators (see e.g. four-loop example in [21]). However, 
for any given algebra the evaluation can be performed easily. It is convenient to present 
the final result for two Wilson lines in the form 



MGEW 


= exp 



^iljl ^*272 

Nr, 



(5.22) 


where E^’s are numeric coefficients. Here, the exponent is the matrix exponent, where 
matrices with indices (fij'i) and { 12 j 2 ) belongs to different matrix spaces. Often one con¬ 
siders the cusp configuration with contracted color indices. Then it is convenient to use 
the following parametrization 



MGEW 


6ij exp 


CX> 

n=l 


(5.23) 


The first few coefficients E for the case of SU{Nc) gauge group are presented in table 1. 
Eor the calculation we have used the ColorMath package [23]. 

In the case of many Wilson lines the general consideration remains the same. Since 
the diagrams contributing to MEK are connected in the absence of Wilson lines they 
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Figure 5. MGEW contributing to MEK for the correlator of four Wilson lines. The diagrams with 
permuted Wilson lines should be added. For light-like Wilson lines these diagrams are zero. 
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Table 2. The coefficients for expansion of correlator of four light-like Wilson lines in MGEW 
approximation (5.25). The variables s, t, u and 9i are defined in (5.26). 


necessarily contain several vertices Vn with n > 1, see the example of two-loop MGEWs 
with four Wilson lines in hg.5, and three-loop MGEWs with three Wilson lines in fig.l (the 
diagrams B,G,E and F). The only non-zero contribution to MEK is given by the one-loop 
diagrams. Therefore, MEK is given by 






MGEW 


E 

k,l=l 

k<l 




(5.24) 


where Ki is given in (5.17), and missed indices should be saturated by Kronecker deltas. 
The expression (5.24) is given for the Wilson lines pointing from the origin. In the case 
of Wilson line incoming to the origin the corresponded generator should be replaced by 
—t^. Thus, the only non-trivial contribution to the correlator of light-like Wilson lines in 
MGEW approximation comes from the defect and can be evaluated algebraically. 

As an illustrative example, we present the correlator of four light-like Wilson lines 
in MGEW approximation in the exponentiated form. The complete result (with eight 
open indices) is rather involved. Therefore, we took a particular convolution of Wilson 
lines, which corresponds to the scattering in t-channel. In MGEW approximation such a 
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correlator is convenient to present in the following form 



MGEW 


(5.25) 


C oo 

n=l 

The expression for coefficients A and B are presented in table 2 up to the three-loop order. 
The variables s, t and u are defined as 



Sn = (^12)”^ + (^34)"^ tn = (^14)"^ + (i;23)"^ + (^24)"^ 

6»i = {V 12 V 13 V 23 Y + (^’l2^^14'^^24)" + {vi3VuV3iY -|- (n23t’24l^34)"■ (5.26) 

The calculation of higher orders in MGEW approximation does not present any diffi¬ 
culties, but also is not interesting. We remind that MGEW is only an approximation, and 
therefore, it can be considered only as an illustration. Nonetheless, the results presented 
in table. 2 show some interesting features. 

One can see that the two-loop expression resembles the one-loop expression. This 
is a known feature of the two-loop soft anomalous dimension, see [24, 25]. Moreover 
the proportionality coefficient is the same as in MGEW approximation for the cusp, see 
E(i 2 ),cusp hr table 1. In fact, it is a consequence of a general statement; in two-loop 
approximation the soft anomalous dimension (for light-like Wilson lines) is equal to the one- 
loop anomalous dimension multiplied by the two-loop cusp anomalous dimension [24, 25]. 
This observation is known as the dipole formula for infrared singularities of scattering 
amplitudes in QCD [26, 27]. The status of the dipole formula at three-loop order and 
higher is currently questionable, e.g. see ref. [28], for recent review see ref. [29]. 

Considering, the expression presented in table 2, we observe that at three-loop order 
MGEW approximation does not support the dipole formula. Since in the considered kine¬ 
matics MGEW approximation is gauge invariant, the expression given in tables 1 and 2 
represent some irreducible part of the complete result. In this way we conclude that the 
violation of the dipole formula at the three-loop order is likely. 


6 Conclusion 

In this article we have presented an approach that allows one to evaluate the correlator 
of a product of Wilson lines directly in the exponentiated form. In contrast to the ex¬ 
isting methods of exponentiation, e.g. [3-7], the presented method does not rely on the 
diagrammatical consideration, and allows to obtain the argument of exponent in the form 
of correlators of known operators. To our best knowledge, it is the first example for such 
a relation in a non-Abelian gauge theory. 

The approach is based on the fundamental property of the perturbative expansion, 
namely, on the exponentiation of the connected part of Eeynman diagrams. In sec. 3.1 
(see also [9]) we show that this property implies the exponentiation of the diagrammatic 
expansion for a wide class of operators; the operators, which can be presented in the form 
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of the exponent (3.3). Although this observation is evident and well-known, it was not 
regularly applied for the consideration of Wilson lines and their correlators. 

In the line of derivation it was natural to split the problem of exponentiation onto two 
independent problems; the exponentiation within the non-Abelian gauge theory, and the 
exponentiation of a matrix object (see sec. 4.2). As the result of such consideration the 
exponentiated expression is split onto two components, those are the matrix exponentiation 
kernel (MEK) and the defect of exponentiation. These two components are very different 
both in the form and in the content. MEK contains all the essential information about the 
exponentiated expression. While the defect represents the difference between the matrix 
and the scalar exponentiation procedure, and it is an algebraic function of MEK. 

In this way, the knowledge of MEK results to the knowledge of the complete exponent. 
The expression for MEK has the form of a sum of correlators of the operators V (4.20). 
The operators V have very peculiar form of nested commutators of gauge helds, and play 
an exceptional role in the non-Abelian exponentiation. Undoubtedly, further investigation 
should reveal exceptional properties of the operators V, and hence, of MEK. Good illustra¬ 
tion of non-trivial properties of MEK is the correlator of light-like Wilson lines in MGEW 
approximation. In this case MEK is given exactly by one-loop diagram (see sec. 5.2). 

The defect of exponentiation is an algebraic function of MEK. In sec. 4.4 we have 
shown that the defect is given by the simple recursive formula (4.30), where the most non¬ 
trivial operation is the procedure of the matrix symmetrization. Therefore, the evaluation 
of the defect can be easily automated. 

In sec. 5, we have performed the comparison of the presented approach with existing 
approaches. We have shown the complete agreement between the presented approach 
and the standard approach to the cusp exponentiation at two-loop order [3, 4] (Here, the 
standard approach is the approach based on the modification of color factors, which can 
be applied only for singlet conhgurations of Wilson lines). We have also demonstrated 
an agreement and the relation of the presented approach to the exponentiated diagrams 
obtained by the replica trick [5, 6]. Remarkably, the splitting of the exponent onto MEK 
and the defect, can be also traced within the replica approach. Thus, we conclude that the 
presented approach is in complete agreement with other methods of the exponentiation. 

Using MGEW approximation we have analyzed the light-like configuration of Wilson 
lines. In this approximation and kinematics the web diagrams can be evaluated with a 
minimal effort, they are given entirely by one-loop MEK and the corresponded defect. 
Remarkably, that MGEW part of the complete expression is gauge invariant that allows 
us to make a judgement about the form of the whole expression. Using this approximation 
we observe the violation of the dipole formula for infrared singularities of parton scattering 
amplitude at the three-loop order. However, we can not strictly conclude the violation of 
the dipole formula due to the possible cancelation of violating terms with similar terms 
from the diagrams do not contributing MGEW approximation. 

The presented procedure of exponentiation allows further generalizations. So, for ex¬ 
ample, the correlator of ordinary product of Wilson lines (i.e. the diagrammatic series with 
exchanges by real particles) can be exponentiated by the same procedure, as it is demon¬ 
strated in sec. 3.3. Therefore, the presented exponentiation procedure can be applied to a 
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wide range of physical tasks, e.g. calculation of the soft anomalous dimension, the thresh¬ 
old resummation, the evaluation of the soft factor for transverse momentum dependent 
factorization, description of diffractive processes. 
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A Operator 14 at low orders 


In this appendix we present the explicit expression for the operators 14 up to the fourth 
order. These expressions are used for the calculation in the sec. 5. 

The operator is defined in (4.1,4.2), while the operator 14)is defined in (4.6). 
Elaborating the matrix structure we obtain the expressions for the first four operators 
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)} 


where A* = 7^(rj)4A“( 7 (Ti)) and Ao = 7 ^(T)taA“( 7 (r)). Here, as everywhere in the 
article, 7 (t) is a parametrization of a Wilson line curve, 7 (t) is a tangent to the curve, 
and ta is the generator of a gauge group. The expressions for 14^4 after the first equality 
symbol have the form resulting directly from eqns. (4.1,4.2), while the expressions after 
the second equality symbol are obtained by rearranging the integrals and renaming the 
integration variables. 
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The expression (A.l) is valid for the Wilson line in any representation. The represen¬ 
tation dependence is concentrated in the generator prefactor of operator V in eqn.(4.2). It 
can be archived by the selection of the proper normalization condition for the generators 
and structure constants. The normalization condition used in (A.l) is 








^ab 


for generators of any representation. 

For practical application one usually consider the straight segments of Wilson lines. 
Let us consider a straight Wilson line from the point to the point z^. In this case the 
contour parametrization reads 


7 (r) = -t-(1 - (A.2) 

and the tangent to the contour is the constant vector 

7 ^(r) = = v^. 


The Feynman rules for the operators Vn in the position space are defined as 
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a,ai... 


"{xi, 


^) = 




-I/“ 

1 7 
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(A.3) 


Therefore, for a straight Wilson segment (A.2) the Feynman rules corresponding to oper¬ 
ators Vi, 2 , 3,4 read 


Va!^al^aUxi,X2) 

K!'alala!ixi,X2,X3) 
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( 6*12 - O21) , 

. 3 
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X (^1234 — ^1243 — ^2341 + 6*2431 “ ^3214 + ^3421 + ^4213 “ ^4321) 
+'^1324 +'F’i432 +'^2134 + '^2431 +'^2341 +'^3142 +'^3412 +'^3214 
+'^4132 + 'P4312 + 'P4213 ) 


where all Xi are located on the path of Wilson line, and 0*,,^ = 6{y < Xi < .. < Xk < z) the 
Heaviside function ordering coordinates along the Wilson line. The symbol Vijk stands for 
the first term in the bracket with permuted color indices and coordinates, i.e. { 010203 } —)• 
{ojOjOfc} and {X 1 X 2 X 3 } ^ {xiXjXk}- 

In the case of half-infinite straight Wilson lines it is convenient to use the same 
parametrization as for the segment (A.2). The only change is that the scalar parame¬ 
ter runs over the infinite range r G (0, 00 ). Such a configuration of Wilson lines can have 
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additional infrared divergences, resulting from the interaction with gluons at infinity. In 
order to regularize these divergences one commonly uses the suppressing exponent with 
infinitesimal argument [22, 30]. Formally such a regularization implies the change of the 
gauge field within the Wilson line 

A^,{vT + y) ^ + (A.5) 

where 5 and e are infinitesimal and 5 > e > 0. Within such a regularization the large- 
distance singularities are regularized for space-, time and light-like Wilson lines. Thus, 
the change of the Feynman rules for a half-inhnite Wilson lines is minimal. Namely, the 
function should be replaced by 9{y < Xi < ... < Xk < oo)exp v‘^ — ie Xi^ . 
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